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Two-dimensional (2D) strongly correlated elec- 
tron systems underlie many of the most impor- 
tant phenomena in contemporary condensed mat- 
ter physics, including the Quantum Hall Effect 
(QHE) [Q, "high Tc" superconductivity [§, and 
possible exotic conducting states in silicon MOS- 
FETs l^. We demonstrate the existence of yet 
another exotic ground state in strongly corre- 
lated, 2D electronic materials: a novel, insulating 
bond-order/charge density wave state (BCDW) 
in the commensurate 1/4-filled band that per- 
sists for all anisotropies within the 2D lattice, 
in contradiction to the non-interacting electron 
prediction of the vanishing of density waves in 
2D for non-l/2-filled bands. The persistence of 
the BCDW in the 2D lattice is a consequence 
of strong electron-electron (e-e) interaction and 
the resultant "confinement," a concept recently 
widely debated 0]. Our results have implica- 
tions for experiments in the organic charge trans- 
fer solids (CTS), where they explain the obser- 
vation of a "mysterious" coexistence of density 
waves pjq], clarify the optical conductivity of the 
"metallic" state [Q, and suggest an approach to 
the observed organic superconductivity. 

We consider a 2D lattice, consisting of coupled chains 
of strongly correlated electrons, described by the quasi- 
2D extended Peierls-Hubbard Hamiltonian 
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varies from to t, the electronic properties vary from 
quasi-lD to quasi-2D, modeling a wide range of materi- 
als. Each site is occupied by a molecular unit, the dis- 
placement of which from equilibrium is described by Uj^ra 



is an intramolecu- 
h.c], where 



In the above, j is a site index, ?7i is a chain index, a 
is spin, and we assume a rectangular lattice W\. As t_L 



(with Aj^„j = {Uj^m - Wj + i^m)), Vj^ra 

lar vibration, and Bj,k,e,m,a = [c]^f^^^Ck,m,a 
c\ I ^\s a Fermion operator. We treat the phonons in the 
adiabatic approximation and are interested in uncondi- 
tional broken symmetry solutions that occur for electron- 
phonon (e-ph) couplings (a, /?) ^ 0+. Although we limit 
our explicit analysis of e-e interactions to on-site repul- 
sion U and intrachain nearest-neighbour interaction V 
only, we will find that additional intersite interactions 
wiU further stabilize the BCDW. 

For the 1/2-filled band, the difference between ID and 
2D is well understood and profound. In the ID limit, 
widely discussed for polyacetylene |g|, the ground state 
for physically relevant values of t/ > 2V > is a bond- 
order wave (BOW), with Uj — (— 1)^mo. In the isotropic 
2D limit, which has been applied to the parent com- 
pounds of the high Tc superconductors, the ground state 
for the same physically relevant parameters is an anti- 
ferromagnet (AFM, i.e., a. 2ki? SDW) @. Thus in the 
1/2-filled band, the dominant broken symmetry depends 
very strongly on dimensionality and (a) there is no co- 
existence between the BOW and the SDW; and (b) the 
strength of the SDW increases monotonically with the 
interchain hopping. 

In contrast to the 1/2-filled band, broken sym- 
metries in non-l/2-filled commensurate correlated 
bands have been investigated chiefly in the quasi-lD 
regime, i^ << t, and have employed primarily Fermi 
surface/" k" -space arguments, which imply that finite t± 
destroys the Fermi surface nesting that characterises the 
ID limit, leading necessarily within the simplest k-space 
interpretation to the restoration of the metallic phase 
|ll| . The nesting concept, however, is provably appli- 
cable only in the limit of zero e-e interactions. For 
strongly interacting electrons, broken spatial symmetries 
are more naturally studied by "configuration space" ar- 
guments [pyi2|]. Using such arguments, we show be- 
low that the ground state of Eq. (1) is a novel BCDW 
state that persists for all t±_ for the strongly interact- 
ing regime U > 4|i|, provided that the nearest-neighbour 
Coulomb interaction V < Vc, where Vc — 2\t\ in the limit 
[/ — » oo and is slightly larger for finite U. For small t±, 



the BCDW state drives a SDW (creating a "BCSDW", 
|^,|j|), but the SDW amphtude, after reaching a maxi- 
mum, vanishes as tj_ is further increased. The persistence 
of the BCDW to the isotropic 2D hmit, its coexistence 
with the SDW for smaU (but nonzero) t±, and the van- 
ishing of the SDW at large t±, are all new results, quite 
unexpected from established behavior of the 1/2-filled 
band. 

We begin with the physical intuition suggesting the 
existence of the 2D BCDW and the peculiar behavior 
of the spins. For the ID 1/4-filled band, Hubbard (l|] 
showed that for large U and in the presence of long-range 
Coulomb interactions, there exist two distinct Wigner 
crystals, represented by the configuration space t = 
"cartoons" as ...1100... and ...1010..., where the num- 
bers denote electron site occupancies (for finite hopping, 
1(0) actually represents site occupancy 0.5 -I- e (0.5 - 
e), with e small but nonzero). Ref. iQ showed that 
long-range Coulomb interactions are not essential for the 
Wigner crystal ..1100.. ; this ground state is also ob- 
tained for the ID limit of Eq. (1), provided V < Vc\ only 
for V > Vc does the the Wigner crystal ....1010.... be- 
come the ground state pq] . Importantly, for nonzero a or 
/3, this ...1100... ground state is a coexisting BOW-CDW 
Jig , i.e., BCDW, as indicated by the individual chains in 
Fig. 1. The bonds in the ID limit are singlets, but as t± 
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FIG. 1. Sketch of the novel BCDW ground state that co- 
exists with the SDW for small t± in the strongly correlated, 
anisotropic 2D 1/4-filled band. The arrows indicate the spin 
directions and their sizes indicate the relative charge and spin 
densities. The hopping integrals used to calculate the energies 
of the distorted lattices correspond to r4 = (see text) and 
are shown above the bonds along the top chain. This variation 
reflects the BOW. Note that the charge ordering corresponds 
to ID Wigner crystals along the longitudinal, transverse, and 
both diagonal directions (see text). 

is gradually increased, the sites acquire spins M], as in- 
dicated in Fig. 1. The spin arrangement along a chain is 
explained by inspection of the ID Wigner crystal cartoon 
...1100....: each '0' is closer to a specific '1', and therefore 
for finite hopping acquires the same sign (though not 
magnitude) of the spin density as the site with larger 
charge density (since it is the same electron that is 
shared between these two sites). Our detailed numeri- 
cal calculations show that a 7r-phase shift between chains 



gives the lowest electronic energy, implying the structure 
shown in Fig. 1 for the t± << \t\ limit. In this Figure, 
the modulations of the hopping integrals are determined 
by Uj,,n = uo[(-l)™'"2Cos(2fcFJ - 7r/4) + r4 cos(4fcFi)], 
where r2 and r^ are the amplitudes of the 2kF and 
Akp distortions of Uj,„i, respectively. Our many-body 
calculations show that this leads to intrachain CDW 
PciJ) = {'Ecr'^]j,<y^j,e,<y) = 0-5+ Pocos(2/CFJ - 37r/4), 
SDW psij) = (c]^ |C-,-^£,| - cli^^Cj^i^i) = ps2Cos{2kFJ - 
7r/4) + ps4(cos(4fcFJ-7r), and BOW (I]^ -^jj+1:™.™:<^) = 
bo + b2COfi{2kpj — 7T/2) + b4COs{Akpj — it) where the phase 
angles of the DW's given correspond to odd m, and their 
amplitudes depend on [/, V, and t± . 

We now discuss the consequences of increasing t± fur- 
ther. Interchain hopping t± leads to partial double occu- 
pancy on a single site (t i) with an energy barrier that, 
while less than the bare U, is a Ueff that increases with 
U. For large enough Ueff, the electrons are thus confined 
to their respective chains, and the BCDW state persists 
up to the isotropic 2D limit for the strongly correlated 
case. A striking feature of the BCDW state is that it 
is a Wigner crystal along the chains (...1100..., periodic- 
ity 2kF), transverse to the chains (...1010..., periodicity 
4ki?), as well as along both diagonal directions (...1100..., 
periodicity 2'kF). The BCDW state is thus a particularly 
robust broken symmetry M. For instance, by enhanc- 
ing the 4kF charge order along the transverse direction, 
Vj_ enhances the stability of the BCDW. Similarly, the 
diagonal ...1100... charge ordering implies that even the 
additions of hopping i^iag and Coulomb repulsion Vdiag 
along the diagonals will not destroy the Wigner crystal 
for realistic parameters (in particular, Vdiag stabilizes the 
BCDW relative to the Wigner crystal that is ...1010... 
along both x and y directions). 

We next describe the evolution of the spin structure. 
From the cartoon in Fig. 1, we see that for the SDW to 
exist it is essential that the 'O's have spin. In the small 
t± case, the sign of the spin on a '0' is necessarily that 
of the nearest intrachain '1'. Note, however, that each 
'0' also has an intercha,in '1' as a neighbour and that 
for a stable SDW the spin densities of the 'I's that are 
neighbours of a specific '0' must be opposite. Therefore, 
with increasing t±, competing effects occur. On the one 
hand, the magnitude of the interchain exchange coupling 
~ t\/Ueff increases. On the other hand, the spin density 
on a site labeled '0' decreases because of the cancelling 
effects of the intra- and mier-chain neighbouring 'I's. We 
thus expect the SDW to vanish at a i]_ that will depend 
on the magnitudes of the bare U and V. 

To confirm these expectations we use exact diago- 
nalization and constrained path quantum Monte Carlo 
(CPMC) ||l^ numerical techniques to calculate for rep- 
resentative finite 2D lattices: (i) the electronic energy 
gained upon bond distortion, AS = E{Q) — E{uj_,n), 
where E{uj^„i) is the electronic energy per site with fixed 



distortion Uj,m along the chains; (ii) site charge densi- 



ties Tij and intrachain bond orders 
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and (iii) the z-z component of the spin-spin correlations, 
for each U, V and t±. A decreasing AE as a func- 
tion of t± signals the destruction of the distortion by 
two-dimensionality, while a constant or increasing AE 
indicates a persistent distortion [B 12 Iq]. In order to 
determine the correct behaviour in the thermodynamic 
limit from finite-size simulations, we choose lattices and 
boundary conditions based on the physical requirement 
that any nonzero t± must destabilise the BCDW for non- 
interacting electrons on that particular finite lattice (see 
supplemental information). We have studied both r4 = 
and r4 = r2 (see supplementary information), estab- 
lishing that while the magnitude of AE does depend on 
r4/r2, its behaviour as a function of t± does not, pro- 
vided only that ^2 7^ 0. For brevity we present here 
the results for r^ = only. Wc note that the magni- 
tudes of AE for the noninteracting and the interacting 
cases are very different: Coulomb interactions reduce the 
AEq = AE \tj_=o considerably. However, this merely in- 
dicates that for a given e-ph interaction the magnitude 
of the distortion is less for correlated electrons than for 
noninteracting electrons. Since, however, the interacting 
single chain is distorted |lq| , and since our interest lies in 
determining the behaviour as a function of t± only, the 
relevant quantity is not the absolute value of AE but the 
normalised energy gained per site upon distortion, i.e., 
AE/AEq. This argument is similar to that made by An- 
derson Q for perturbative treatments of coupled chains, 
viz., the sequence in which the intrachain Coulomb inter- 
actions and interchain hopping are included is important, 
and a correct physical picture is obtained only by first in- 
cluding the Coulomb interactions. 

In Fig. 2(a) we compare the behaviour of AE/AEq 
for the non-interacting and interacting {U — 6|i|, V — 
\t\) cases for three different lattices satisfying our bound- 
ary condition constraints. Finite 4n-electron ID periodic 
rings have a strong tendency to have total spin S = 1 for 
nonzero U. In the present case, the 8-site periodic ring is 

5 = 1 in both the undistortcd and distorted states, but 
the 16-site ring is S = 1 in the undistorted state and S 
= when distorted. The AEq for the 8x2 and 8 x 

6 then correspond to the energy gained upon distortion 
by the ID S = 1 state. This is appropriate, since the 
ground states of the 8x2 lattice are S = for the small- 
est nonzero ij_, and as seen in Fig. 2(a), the exact AE 
for this case converges smoothly to the AEq at ij^ ^ 
(the exact diagonalizations here included calculations for 
i_L as small as 0.01). For the 16 x 6 lattice, the one- 
electron orbital occupancy at [/ = is degenerate even 
for the distorted state in the region t\_ > 0.6 indicating 
that the distorted state ground state is S = 1 for large 
t± and nonzero U. The calculated AE for the 16x6 are 
therefore compared to the triplet A_Eo in the region t±_ 
> 0.6. 



8x2 U=V=0 
8x6 U=V=0 
16x6 U=V=0 
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FIG. 2. AE/AEo for the 8 x 2, 8 x 6, and 16 x 6 lattices 
as a function oit±, for U = 6t,V — t. (b) Site charge densities 
on the sites along a chain of the bond-distorted 8x6 lattice, 
for different t± and the same U and V as in (a). 



As seen in Fig. 2(a), while for the non- interacting cases 
the AE/AEq decreases rapidly with t±, for the interact- 
ing cases the AE/AEq either remains unchanged or is 
enhanced by t±. In Fig. 2(b) we show the charge den- 
sities on the sites along a chain for the bond-distorted 8 
X 6 lattice, as a function of t±, for the interacting case 
only [as expected from the AE plot, the CDW ampli- 
tude decreases rapidly for the noninteracting case]. The 
CDW behaviour is the same for the 8x2 and the 16 x 
6 lattices. Both Figs. 2(a) and (b) indicate the enhanced 
stability of the BCDW in the interacting 2D case. The 
same conclusion is also reached from calculations of the 
intrachain bond orders. For the 8x2 lattice, we have re- 
peated the exact calculations also with nonzero Vj_ , with 
the same conclusions (see supplementary information). 
These results provide quantitative proof of our qualita- 
tive arguments establishing that the BCDW is a robust 
broken symmetry state for the interacting 2D 1/4- filled 
band. 

In Fig. 3 we show the interchain spin-spin correlations 
for the distorted 8x6 lattice for several values of t± . The 
spin-spin correlations indicate a SDW that is enhanced 
between t± =0.1 and 0.4 but then vanishes at t± ~ 0.6. 
We observe this same behaviour of the SDW on 8 x 2 
and the 16x6 lattices. In all cases, the SDW amplitude 
initially increases, exhibits a maximum, and then van- 
ishes at a t^ which decreases with the size of the system. 
The initial increase of the SDW amplitude indicates that 



f]^ is nonzero, a conclusion that is also in agreement with 
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FIG. 3. The z-z spin correlations between sites 1 and 2 on 
the first chain of the 8x6 lattice and sites 1 - 8 on the second 
chain, with U = 6|i|, V — \t\. AFM correlations increase with 
t± up to t± — 0.4 but then vanish at tj_ ~ 0.6, even though 
the BCDW continues to persist for all t± (see Fig. 2). 

the experimental observation of the BCSDW state in the 
weakly 2D organic CTS ||,^. Based on the calculations 
for 16 X 6 lattice, we estimate 0.1 < f]^ < 0.3. Theo- 
retically, the 2D BCDW introduces a new mechanism for 
broken symmetry in strongly correlated lattice systems. 
The persistence of the BCDW state at large t± is related 
to the Wigner crystal-like charge arrangement along all 
possible directions. Indeed, our 1/4-filled BCDW state is 
remarkably similar to the "paired electron crystal" state 
found in the continuum electron gas by Moulopoulos and 
Ashcroft m 

Experimentally, organic CTS, which span the range 
t± < O.lt in (TMTTF)2X to t± ~ i in certain (BEDT- 
TTF)2X, provide a critical testing ground for our re- 
sults. First, the 1/4-filled band and the resulting co- 
existence between the BCDW and the SDW at small 
t± provide a natural explanation [fsUTj] for the other- 
wise "mysterious" coexisting CDW/SDW state that has 
been observed in (TMTSF)2PF6 and (TMTTF)2Br |] 
and likely also in a-(BEDT-TTF)2MHg(SCN)4 |]. Sec- 
ond, recent data M on frequency-dependent optical con- 
ductivity in the high temperature "metallic" phase of 
(TMTSF)2C104 and (TMTSF)2PF6 exhibit both zero- 
and finite-energy modes, with the zero-energy mode hav- 
ing a spectral weight of only about 1 %. Assuming 



that the (TMTSF)2X are weakly incommensurate 0, 
our BCDW state provides a natural explanation, with 
the zero-energy mode ascribed to the conductivity of the 
incommensurate defects and the gap mode to the excita- 
tion across the BCDW pseudogap that persists at high 
temperature due to fluctuations associated with low di- 
mensionality. 

What might be the relation of our new BCDW state 
to organic superconductivity, the mechanism for which 
remains unclear despite two decades of research? Experi- 
mentally, in the (TMTSF)2X, the BCSDW state for small 
t± gives way to superconductivity under pressure p9|. 
Insulator-superconductor transition has also been noted 
in the BEDT class of materials, where superconductiv- 
ity is often obtained by slight modifications of the anion 
(for instance, a-(BEDT-TTF)2MHg(SCN)4 are BCSDW 
for M = K, Rb and Tl, while the M = NII4 compound 
is superconducting; k-(BEDT-TTF)2Cu[N(CN)2]X is an 
antiferromagnetic semiconductor for X = CI, but X — Br 
is superconducting) ^^. Both pressure and larger an- 
ions are believed to enhance the t± . One-electron theory 
might suggest that the appearance of superconductivity 
is related to the complete destruction of the background 
DW. Thus the persistence of the BCDW state, even in 
a region where the SDW has vanished, suggests a dif- 
ferent scenario: namely, interconnected or coexisting di- 
agonal and off-diagonal long-range order, perhaps in the 
presence of a few commensurability defects. Such a sce- 
nario has been discussed by Moulopoulos and Ashcroft 
in their work on the paired electron crystals [H, and 
goes back to 1970 uM. Similarly, the possibility of an 
insulator-superconductor transition in the dilute 2D elec- 
tron gas, where the insulator is once again a "paired" 
Wigner crystal, has recently been discussed p^. Al- 
though the "paired" Wigner crystal in the 2D electron 
gas remains to be proved and is presently controversial 
0, it is intriguing that the BCDW state found here is ob- 
tained naturally within the standard microscopic lattice 
Hamiltonian for organic CTS. Several features of our re- 
sults support the above speculation about the insulator- 
superconductor transition. First, the vanishing of the 
SDW and the persistence of the BCDW, taken together, 
suggest the possibility of singlet coupling between the 
sites labeled '1' in Fig. 1 in this region, leading to a spin 
gap. Second, within many existing models of supercon- 
ducting pairing involving correlated electrons, the inter- 
actions that bind two particles often also lead to phase 
segregation. In contrast, any binding of commensurabil- 
ity defects here would be a consequence of the "quasi- 
chemical bonding" that occurs between the "occupied" 
sites, and these pairs are already separated. Finally, our 
preliminary calculations in the region t± > i^. in small 
lattices show signatures of superconducting pairing cor- 
relations (see supplementary information), but these cal- 
culations will have to be extended to larger lattices before 
firm conclusions can be drawn. 
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